Envelope quasi-solitons in dissipative systems with cross-diffusion 



V. N. Biktashev 

Department of Mathematical Sciences, University of Liverpool, Liverpool L69 7ZL, UK 



o 

< 

<N 



(N 
> 

m 
o 
o 
cn 

o 



M. A. Tsyganov 

Institute of Theoretical and Experimental Biophysics, Pushchino, Moscow Region, 142290, Russia 

(Dated: January 21, 2013) 

We consider two- component nonlinear dissipative spatially extended systems of reaction-cross- 
diffusion type. Previously, such systems were shown to support "quasi-soliton" pulses, which have 
fixed stable structure but can reflect from boundaries and penetrate each other. Presently we 
demonstrate a different type of quasi-solitons, with a phenomenology resembling that of the envelope 
solitons in Nonlinear Schroedinger equation: spatiotemporal oscillations with a smooth envelope, 
with the velocity of the oscillations different from the velocity of the envelope. 

PACS numbers: 87.10. +e, 02.90.+p 



INTRODUCTION 

Dissipative structures, i.e. patterns in spatially ex- 
tended systems away from equilibrium have been inten- 
sively studied for many decades. A very comprehensive 
review can be found in Cross and Hohenberg [1]; re- 
sults obtained since then would probably require an even 
more extensive review. A very popular class of math- 
ematical models is the reaction-diffusion systems with 
diagonal diffusion matrices. There have been numerous 
indications that non-diagonal elements in diffusion ma- 
trices, i.e. cross-diffusion, can lead to new nontrivial ef- 
fects not observed in classical reaction-diffusion systems, 
e.g. quasi-solitons in systems with excitable reaction part 
[2, 3]. The defining features of the quasi-solitons was 
their ability to penetrate each other, which makes them 
akin to the true solitons in the conservative systems. 
However the question remained whether this similarity 
is a reflection of common mechanisms, or is entirely su- 
perficial and incidental. Here we report an observation 
which makes the similarity even more striking. Namely, 
the previously reported quasi-solitons propagated while 
retaining fixed shape profile, i.e. were constant solu- 
tions in a CO- moving frame of reference; the exceptions 
were the "ageing" quasi-solitons reported in [3] which re- 
tained their front and tail structures but changed their 
overall length. Here we report "envelope quasi-solitons" 
(EQS), which share some phenomenology with envelope 
solitons in the Nonlinear Schroedinger Equation [1, 4]. 
Namely, they have the form of spatiotemporal oscilla- 
tions ("wavelets") with a smooth envelope, and the ve- 
locity of the individual wavelets (the phase velocity) is 
different from the velocity of the envelope (the group ve- 
locity). This may be a serious evidence for some deep 
relationship between these phenomena from dissipative 
and conservative realms. 

Our observations are made in two two-component 
models, supplemented with cross-diffusion, rather than 
self-diffusion terms; such terms may appear say in me- 



chanical [5], chemical [6], or ecological [7, p. 11] applica- 
tions: 
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We consider the FitzHugh-Nagumo (FHN) kinetics taken 
in the form 



/ = u(u — a)(l — u) — k\v^ 



(2) 



as an archetypal excitable model, with an arbitrarily 
fixed value of parameter k\ = 10, and varied values of 
parameters a and e. As a specific example of a real-life 
system, we also consider the Lengyel-Epstein [8] (LE) 
model of chlorite-iodide-malonic acid-starch autocatalitic 
reaction system. 
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for fixed A = 6.3 and B = 0.055. We simulated (1) on 
an interval x G [0,L], L < oo with Neumann boundary 
conditions for both u and v [9]. 

Fig. 1 illustrates development and subsequent propaga- 
tion of an envelope quasi-soliton (EQS) solution in (1,2). 
The profiles are presented in a co-moving frame of ref- 
erence, with x-coordinate measured with respect to the 
"centre of mass" Xc of the quasi-soliton [9]. Simulations 
with different initial conditions show that as long as the 
initial perturbation is above a threshold, the amplitude 
and overall shape of the quasi-soliton does not depend 
on its details. An important feature, evident from com- 
paring panels (d) and (e), is that whereas the overall 
shape (the envelope) of the quasi-soliton and the wave- 
length of the high-frequency oscillations (the wavelets) 
within that envelope remain unchanged, the phase of the 
the wavelets relative to the envelope position changes, so 
the phase velocity (of the wavelets) is different from the 
group velocity (of the envelope). 

This feature can also be seen in fig. 2(b). The thin 
stripes in the density plot represent individual wavelets. 
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FIG. 1: (color online) Quasi- soliton profiles at the indicated 
moments of time, in a co-moving frame of reference, upper 
x-axes, with the reconstructed original spatial coordinates 
shown on lower x-axes. Parameters a = 0.12, s = 0.01, 
L = oo, solution propagates rightwards. (a-c) Develop- 
ment of a quasi-soliton. (d,e) Propagation of a quasi-soliton, 
with unchanged envelope and shifting phase of high-frequency 
wavelets within the envelope. Note that wavelets in (d) and 
(e) are in antiphase: the v profile at x = Xc is near a local 
minimum in (d) and a local maximum in (e). 



and the broader band, consisting of these stripes, repre- 
sents the envelope. The slope of the stripes is the inverse 
of the phase velocity, and the slope of the band is the 
inverse of the group velocity. The stripes are not parallel 
to the band, because the group and the phase veloci- 
ties differ. This figure also illustrates another important 
phenomenon: the reflection of the quasi-soliton from the 



boundary. 

Panels (a) and (c) in fig. 2 illustrate two different sorts 
of solutions which are observed at higher and lower values 
of parameter a, which also reflect from the boundary. 

In panel (a), we still see individual wavelet stripes that 
are not parallel to the envelope bands, but there are two 
bands in the incident wave. Note that the reflected wave 
only has one band; however if that reflected band is al- 
lowed to propagate for a sufficiently long time, it will 
spawn its twin band behind it. This is a "multiplying" 
EQS. We do not go further into properties of these solu- 
tions, reserving that for a future study. 

In panel (c) there is only one dominant stripe and many 
weaker stripes, all of which are parallel to the band. This 
solution has phenomenological features similar to quasi- 
solitons described previously, e.g. in [2], namely, the wave 
retains constant shape as it propagates, and reflects from 
a boundary. 

Panel (d) shows a quasi-soliton reflecting from the 
boundary, in the other model (1,3). 

Fig. 3 gives an overview of the parametric area of the 
EQS solutions in (1,2) and its neighbours. In panel (b), 
the parameter sets at which EQS solutions like the one 
shown in fig. 2(b) have been observed, are designated by 
red solid circles. This area is surrounded: 

• at higher and lower values of £, by solutions which 
have similar shape to those shown in fig. 1 and 
fig. 2(b), but do not reflect from boundaries ('an- 
nihilating', blue crosses); 

• at lower values of a, by multiplying EQSs, one of 
which is illustrated in fig. 2(a) ('multiplying', green 
stars); 

• at higher values of a, by constant shape quasi- 
solitons, such as the one shown in fig. 2(c) ('con- 
stant', magenta triangles). 

The area of existence of all these solutions in the (a, s) 
parametric plane is bounded from above and from the 
right, and beyond it our initial conditions did not produce 
any stably propagating solutions ('decaying', black open 
circles). Panel (a) in this figure illustrates the variabil- 
ity of the shape of EQS within their parametric domain. 
Most important conclusion from fig. 3 is that the EQSs 
are not a unique feature of a special set of parameters 
but are observed in a rather broad parametric area. 

The oscillatory character of the fronts of cross-diffusion 
waves, described in numerical simulations [2, 3] and anal- 
ysed theoretically in [2, 3, 10], was for waves of station- 
ary shape. Although the proper theory of the EQSs is 
beyond the scope of this Letter, the analysis of their 
non-stationary front structure is easily achieved via lin- 
earization of (1). The resting states in both FHN (2) 
and LE (3) kinetics are stable foci which already shows 
propensity to oscillations. Considering in more detail the 
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FIG. 2: Density plots of the quasi-solitons reflecting from a boundary [9]. (a-c) FHN kinetics, u- = —0.3, = 1, e = 0.01, 
and a is varied as shown under the panels, (a) A double EQS becomes a single EQS upon the reflection. Some time after that, 
it will grow its twin become a double quasi-soliton again, (b) A single EQS: this is the same case as shown in fig. 1. (c) A 
"classical" quasi-soliton which retains its shape as it propagates, (d) An EQS in the LE kinetics, L = 300, u- = 0.8, = 3.3. 
In all panels, the origin of the t-axis is shifted to an arbitrarily chosen moment shortly before the impact event. 
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FIG. 3: (color online) (a) The number of wavelets in an EQS 
as a function of a at fixed s = 0.01. (b) Areas of different sorts 
of solutions in the parametric plane (a,s). The black dashed 
line corresponds to the parametric cross-section shown in (a). 



FHN kinetics, for small u^v^ the solution has the form 



(4) 
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FIG. 4: (color online) Profiles of an EQS wavefront and its 
fitting by (4) at selected moments of time. Parameters are 
£ — 0.01, a — 0.12, L — oo. The origin of the x-axis is chosen 
arbitrarily. 



where 



A(A,i/)v = 0, VT^O, 
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det A = 0, 

/ic — iuj^ V = 
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(5) 
ik. 



Fitting of the ^-component of a solution shown on fig. 1 
to (4) gives c ^ 4.07698, k ^ 1.71532, /i ^ 0.182305 and 
u ^ 6.15190, which satisfies (5) to 3 s.f. [9] The quality 
of the fitting is illustrated in fig. 4. Note that the phase 
velocity of the wavelets here is Cph = u/k ^ 3.59, smaller 
than the group velocity, c ^ 4.08, which agrees with the 
fact that the slope of the individual stripes in fig. 2(b) 
(which is the inverse of the phase velocity Cph) is steeper 
than the slope of the band (which is the inverse of the 
group velocity c). 

The shape of the profiles in fig. 1 is reminiscent of local- 
ized states in the generalized Swift-Hohenberg equation 
with "snakes and ladders" bifurcation diagrams [11]. The 
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essential difference of our solutions is that they move and 
do not preserve their shape, so cannot be immediately 
studied by ODE bifurcation techniques. 

The defining features of the EQSs described above are 
similar to envelope solitons of the Nonlinear Schrodinger 
Equation. The version of this equation known as 
'NLS+' [4], can be written in the form 







for a complex field which presents a reaction-cross- 
diffusion system for two real fields u and v via w = u — iv 
of the form (1) with 



(6) 



System (1,6) has soliton solutions in the form of (fast) 
harmonic waves with a unimodal (sech-shaped) envelope, 
and the propagation velocity of the envelope (the group 
velocity) is different from the propagation velocity of the 
wavelets (the phase velocity). Hence one might think of 
possible interpretation of the EQSs in (1,2) or (1,3) as a 
result of a non-conservative perturbation of the envelope 
solitons in (6), which would select particular values of 
the otherwise arbitrary amplitude and speed of the soli- 
ton and modify its shape. This interpretation, however, 
does not seem to work, and our attempts to connect the 
solutions in (1,6) and (1,2) via a one-parametric family 
of systems have been unsuccessful, as the EQS solutions 
disappeared during parameter continuation. The appar- 
ent reason is that the sense of rotation of solutions of 
(6) in the (tx, v) is clockwise whereas in (2) and (3) it 
is counterclockwise, and the variant of (6) with counter- 
clockwise rotation, the 'NLS-' equation, does not have 
envelope soliton solutions. 

Another comparison can be made with "wave packets" 
reported by Vanag and Epstein in microemulsion BZ re- 
action, and corresponding mathematical models, associ- 
ated with finite-wavelength instability of an equilibrium 
in a reaction-diffusion system with unequal self-diffusion 
coefficients [12]. They considered two distinct types of so- 
lutions: small- and large- amplitude wave packets (SAWP 
and LAWP), both capable of refiection from boundaries. 
SAWP are observed in the nearly-linear regime, they have 
the phase speed (of the wavelets) different from the group 
speed (of the envelope, or the packet). However, being 
near to a linear instability and having no stabilizing ef- 
fect of the dispersion as in NLS+, the packets slowly grow 
both in amplitude and in width, i.e. they are not quasi- 
solitons. The LAWP, on the contrary, have fixed ampli- 
tude and width, but their phase and group velocities co- 
incide, so they retain constant shape. They are therefore 
phenomenologically similar to the quasi-solitons reported 
in excitable systems with cross-diffusion [2]. Note that 



adiabatic elimination of a fast component in a reaction- 
diffusion system with very different self-diffusion coeffi- 
cients is one of the ways in which cross-diffusion terms 
may appear [6, 13], so this analogy deserves further in- 
vestigation. 

To conclude, the solutions we have reported resemble 
NLS+ envelope solitons by their morphology and by their 
ability to refiect from boundaries, however they are dif- 
ferent in that amplitudes and speeds of NLS solitons de- 
pend on initial conditions, while in (1) they are fixed by 
parameters of the models. The reported solution are sim- 
ilar to quasi-solitons reported earlier in that they share 
the fixed amplitude and refiection properties, but differ- 
ent in that they do not preserve constant shape as their 
phase velocities are different from their group velocities. 
Hence we believe this is a new nonlinear phenomenon not 
seen before. The mechanisms behind the key properties 
of this new type of solutions require further investiga- 
tion, however it is already clear that this is not simply a 
non- conservative perturbation of NLS. 
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